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Abstract 

Applying the technique of characteristic functions developped for one- 
dimensional regular surfaces (curves) with compact support, we obtain the 
distribution of hitting probabilities for a wide class of finite membranes 
on square lattice. Then we generalize it to multi-dimensional finite mem- 
branes on hypercubic lattice. Basing on these distributions, we explicitly 
construct the Brownian self-transport operator which governs the Lapla- 
cian transfer. 

In order to verify the accuracy of the distribution of hitting probabil- 
ities, numerical analysis is carried out for some particular membranes. 



Introduction 

We are engaged in the physical problems that can be mathematically modelized 
by the Laplacian transfer problem jj. These are : the diffusion through semi- 
permeable membranes and the electrode problem , the heterogeneous catalysis 
||, the NMR in porous environment Q, etc. 

Consider the diffusion of particles from the source (usually supposed planar) 
to the semi-permeable membrane with rather complex geometry, see Fig. ^a. In 
the steady-state regime, the concentration of particles C obeys the Laplacian 
equation, 

AC* = 0. (1) 

The flux of particles in the bulk obeys Fick's law, <I> = — DVC, where D is the 
diffusion coefficient. The flux across the surface is given by <&„ = — WC, where 
W is the permeability of the membrane (the probability per unit time, surface, 
and concentration for a particle to cross the membrane). Equating these two 
fluxes, we obtain the mixed boundary condition, 

often called also Fourier or Robin boundary condition. The physical parameter 
A = D/W plays an important role in this treatment. 

1 E-mail: Denis.Grebenkov@polytechnique.fr 
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Figure 1: The problem of the Laplacian transfer across a resistive irregular 
interface (a membrane or an electrode). 



In the electrode problem (see Fig. |l|b), there are a planar and an irregular 
electrodes (with an interface resistivity r), placed in an electrolyte with a resis- 
tivity p. The electric potential V obeys the Laplace equation Ay = in the 
bulk of the electrolyte. The boundary condition (||) is obtained by equating the 
current from the electrolyte, — VV/p, to the current — V/r crossing the electrode 
surface, with A = r/p. 

For the heterogeneous catalysis, a catalyst of a complex geometry (with 
the reactivity K) is placed in a solution. When a molecule A diffusing in the 
solution, reaches the catalytic surface, it transforms into A* with reaction rate 
K . As above, the concentration of molecules A obeys the Laplace equation 
(0). The boundary condition (||) with A = D/K is obtained due to the mass 
conservation. 

The same arguments allow to describe the NMR in the porous environment 
by the Laplace equation under the mixed boundary condition. So, all these 
different phenomena are described by the same mathematical formalism (see 
§, § for details). 

The Laplacian transfer problem is significantly more complex than the cor- 
responding Dirichlct or Neumann problem. In order to find the solution, Filoche 
and Sapoval || proposed the following program : 

- to choose the appropriate discretization; 

- to solve the discrete problem; 

- to take the continuous limit. 

For these purposes, they introduced the Brownian self-transport operator Q. 
This operator controls the properties of the Laplacian transfer and depends only 
on the geometry of membrane's boundary. It can be easily defined for a given 
discrete membrane: Qk, n is the probability of first contact with the membrane 
on the n-th site if started from the A;-th site (without touching other sites!). 
This is exactly the problem of the hitting probabilities on the lattice. So, if 
one can calculate Q for any "reasonable" membrane, the Laplacian transfer 
problem would be solved (except the difficulties related with the continuous 
limit). In this paper, we construct the Brownian self-transport operator for the 
membranes with a rather general geometry. 

Now we should specify what the term "reasonable membrane" means. For 
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Figure 2: Corners: (a) internal corner; (b) corner point is carried from bound- 
ary to the bulk of membrane; (c) external corner; (d) corner point is removed, 
two links from boundary to corner point are also removed, two other links are 
connected to each other. 



the moment, we confine ourselves to the two-dimensional case (square lattice) 
in order to clarify the representation of our treatment. When two-dimensional 
results are obtained, we shall explain how they can be generalized for the multi- 
dimensional case. On a square lattice, the Laplacian equation (|l|) becomes 

^X,y ~ C X + l^y — C X -l_y — Cx^y+l — C X ^y-\ = 0, (3) 

i.e., diffusion can be modelized by simple random walks. We suppose that the 
membranes are quite regular. It means that one can choose such discretization 
that the number of singular points is negligible by comparison with the number 
of regular points (contrariwise to the case of fractal membranes) . As we decided 
to work on the square lattice, we should represent the membrane as a sequence 
of horizontal and vertical segments. The last difficulty concerns with corners. 
The problem is that a corner point is connected with two lattice points whereas 
a regular point is connected with only one lattice point. We would like to be 
indemnified from these ambiguities^} To avoid the corner points, we use the 
following rules^j: 

- for the internal corner (see Fig. ^|a, |2)b) , we carry the corner point from 
the boundary to the bulk of membrane. 

- for the external corner (see Fig. ^d), we completely remove the corner 
point from the lattice. The links from the boundary to these points are also 
removed. Two other links are closed to each other. 

This operation preserves the connectivity of the lattice and eliminates the 
singularities (corner points). 

Section |l| introduces definitions and restrictions which allow us to reach final 
results. Section || is devoted to the technique developped in our previous paper 
[0 for one- dimensional regular surfaces with compact support. Note that the 
actual treatment is partially based on the previous results. The distribution 
of hitting probabilities for a two-dimensional finite membrane is obtained in 
Section |^. The generalization of these results for multi-dimensional case is 
described in Section |[ In Section |[ we explicitly construct the Brownian self- 
transport operator. Some important generalizations are given in Section ^|. 

2 The other reason to avoid the corner points is related with numerical simulations. Indeed, 
to verify our analytical results we should solve the problem of Laplacian transfer numerically. 
And here it is not clear how to treat the corner poirits. 

3 These rules were described, in particular, in [El. 
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Section |j describes some numerical results in order to verify the method. In the 
last section we make conclusions. 



1 Definitions 

Consider a square lattice Z 2 on a plane. A path is a sequence (finite or infinite) 
of points {A n } c Z 2 such that Vn \A n ~ A n -\ \ = 1, where | A — B\ denotes the 
distance between two points A and B. 

The sequence of points S — {{x m , y m ) G Z 2 : m £ Z} is called membrane's 
boundary if it obeys the following conditions : 

1. Bijection: Index m enumerates the points {x mi y m ), i.e., there is one-to- 
one correspondance (bijection) between the set of integer numbers Z and the 
set of boundary points S. 

2. Boundary: S separates the lattice Z 2 in two disjoint^ sets - £ and M.q - 
external and internal points. 

3. Accessibility: Any point of S is accessible from all external and internal 
points, i.e., for any (x,y) ^ S and any {x m ,y m ) € S there exists a path {A n } 
such that Ai = (x, y) and {A n } n5 = {x m ,y m ). 

4. Compactness: There is only finite number M of boundary points (in S) 
which do not lie on the horizontal axis, i.e., the non-plane part of the membrane 
has a finite size. 

Set M — M.q U S is called finite membrane. 

Once chosen, points (x m , y m ) are assumed fixed in all following calculations. 
We always use notations x m (or x n ) and y m (or y n ) for abscissae and ordinates 
of boundary points. 

In order to simplify expressions and to avoid possible ambiguities, we intro- 
duce the following conventions: 

- the bulk of membrane is placed in the lower half plane (except, possibly, 
a finite number of points); 

- non-plane boundary points are enumerated by index m — 1..M, i.e., 

Vme[l..M] y m ^0, Vm£[l..M] y m = 0. 

Let us discuss the definition of the finite membrane. The boundary and ac- 
cessibility conditions provide that boundary points take their places sufficiently 
close to each other, i.e., Vm 3n \ (x m ,y m ) — (x n ,y n )\ < \/2. The accessibility 
condition prohibits the existence of corner points (see Fig. ||a, Hja) . Moreover, 
this condition forbids also the "diagonal" points (see Fig. |c). This means that 
membrane's boundary is composed of horizontal and vertical segments just as 
required. So, we can conclude that these two conditions provide all the necessary 
properties which were described in the Introduction. 

The compactness condition is essential to achieve our goal. This condition 
means that the membrane is composed by three parts: two plane "tails" on 
the same height, and one intermediate (non-trivial) part that can be very com- 
plex but finite. The importance of the compactness condition was throughly 
discussed in the previous paper, and we do not repeat it here. Note that in 

4 I.e., there is no path from one set to another which does not pass through a point in S. 
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Figure 3: The accessibility condition forbids the corner points (Fig. ||a, ^b) 
and "diagonal" points (Fig. ||c). 



practice, this condition is not too restrictive because usually one considers the 
finite number of sites. 

According to the bijection condition, we can define the function 3{x, y) 
which gives the index of the boundary points (x,y), i.e., J^(x m ,y m ) = in. Note 
that this function is defined only on the boundary points: if (x, y) (fc S, J{x, y) 
has no value. We shall use the following convention: 

if (x, y) ^ iS, the object containing 3{x, y) is equal to 0. 

For example, if we write e l ^( x > v > 6 , we always mean e % ^( x,v ' e Xs{ x -,y)i where 



Xs(x,y) 



1, if (x,y) G S, 
0, if (x,y)$S. 



Note that it is just a useful convention to simplify the expressions. 

We call all the points {(x, y) G Z 2 : y = n} the n th level. We say that the 
membrane's boundary lies between (—N*) th and N th levels if 

N = max{y TO }, N* = - min{y m }. 

The external point A is called near-boundary point, if there exists m such that 
\A — (x m , ym)\ = 1, i.e., A lies "near" the membrane's bourdary. The functions 
defined on these points, are called near-boundary functions (see below). We 
enumerate these points (and functions) by the same index m as for boundary 
points. Note that usually we suppose m € [1,M], i.e., we take into account 
near-boundary points which lie near the non-plane part of the membrane. 

The external points (a;, 0) ^ M. are called ground points. The functions 
defined on these points, are called ground functions. Let us enumerate ground 
points (x g ,0) by index g G [1,G] using bar notation for their abscissae. Note 
that there exist membranes without ground points (see Fig. ||a for example). 

We introduce the outer normal {Sx m ,Sy m ) for each site of the membrane's 
boundary. Due to the accessibility condition, we have no corner points, therefore 
this vector is correctly defined. Note that m-th near-boundary point is just 
(x m + 8x mi y m + Sy m ), i-e., the outer normal is a vector directed from the 
boundary point to the corresponding near-boundary point. 

Let P x , y (n) be the hitting probability, i.e., the probability of the first contact 
with the membrane on point (x n , y n ) if started from an external point (x, y) G £ 
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without touching other boundary points. Their characteristic functions (j)x,y{9) 
are 

oo 

771= — OO 

The inverse Fourier transform allows to obtain P x y (m), 

TT 

P x>y (m) = J ^-e- im6 4> x , y {e). 



2 Technique developped in [|7j] 

Here we briefly present the technique which was developped in [Q] to obtain the 
distribution of hitting probabilities for a regular surface with compact support. 
All the details can be found in this paper. 

The discrete Laplacian equation (^) in terms of characteristic functions is 
simply 



(0) + <Px-l,y <Px,y+l (0) 



(4) 



Using the convolution properties of hitting probabilities and their normalization, 
one obtains the distribution of hitting probabilities for a planar surface, 

7T 

/df) 
_ e ^-x m )e v \y\ {e) = H y_ Xm: (5) 

— 7T 

where function if (6) is 



tp(6) = 2 - cos - v/(2- cos6») 2 - 1. (6) 

We wrote </?' y '(#) given that for the lower half plane (y < 0) one has exactly the 
same result due to reflection symmetry. 

For the non-planar membrane, relations (^J) remain valid for external points, 
but there are also membrane 's points where it is not true. Nevertheless, we 
preserve this form by adding the correction term which equals to on external 
points. 

In order to simplify manipulation with characteristic functions, we introduce 
vectors ^^(9) containing all 4> Xl y(9) on the y-th level with \x\ < L (at the end 
of calculations L — * oo). We can write (||) as 

A&y\9) = ^ y -^(e) + $ fe+1) (6») + A$ fe) (<9), 

where A is tridiagonal matrix with elements: Aij = 4, Ai^+i = Ai+xj = — 1 (wc 
also introduce an insignificant modification A—l l = Al,—l = — 1 to have the 
cyclic structure of A), A$' a ^ is the correction for membrane's points. Matrix A 
has eigenvalues 

\ h = 4-2cos9 h , with 9 h = 2irh/{2L + 1), h e {—L, —,L} 
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and eigenvectors Vh whose Z-th component is simply e lWh . We decompose $^ (9) 
and A$M (0) on the base of 14, 

L L 

Their coefficients c y and Ac y obey the recurrence relations 

X h c y {9, 9 h ) = cy-x(0, 9 h ) + c y+1 (9, 9 h ) + Ac y {9, 9 h ). (7) 

So, the problem now is to find c y (9,9h)- It can be solved in two steps. First, 
using (0) with certain conditions, we can express c y in terms of Ac y , and 
ip(9). Second, we have to find Ac y . The first step was made in 0: 

- taking sufficiently large N u and Ni , we impose two conditions to close the 
recurrence relations (0), 

c Nu+1 (o,e h )mtp(9)c Nu (e,e h ), c- Nl (B,e h ) = o. (8) 

Note that the first condition is an approximate relation; 

- we find the explicit solution of (0), i.e., we express c y in terms of cq and 
{Ac y '} (see Conclusions for more details); 

- we take the limit N u — > oo and Ni — ► oo to obtain 

N 

c y (9,9 h ) = <pM(9 h )co(9,9 h )+ [l y y \0h)+l { - y v \e h )]Ac y ,{9,9 h ), (9) 

y> = -N* 

where 

min{y ,y' } 

7?(^)= E M0h)] 2j - 1+ly - y ' 1 (10) 

3=1 

(here we use the convention that Y^j=a 1S e 1 ua l to if < a, i.e., 7^, is equal 
to if y < or y' < 0). Note that using only two simple identities, 

<p 2 (9 h )-\ h(p (9 h )+i = 0, i y ^ 1 \e h )-\ hl ( y !\e h )+ ly v r 1) {9 h ) = -§ y , yl (y > 0), 

(11) 

one can verify directly, without intermediate steps of , that (^) is the solution 
of recurrence relations (^). Now we write <j>x,y{9) as 



L ix g h r N 

MQ = E if— I <P W (0 h )co(9,9 h )+ (-y?(6 h )+l ( -- y) (e h ))Acy,(9,9 h ) 

h=-L L y ' = _JV* 

(12) 

So, the first step of our program is completed. Note that all the above 
results are exactly the same as for a regular surface with compact support (by 
this reason we did not explain all details of calculation). Now we are ready to 
pass to the second step : to find Ac y >(9,9h) and co(0, 0fi). These coefficients 
significantly depend on the membrane's geometry. Thus we cannot use the 
expressions of [Q, and we should recalculate Ac y >(6,6h) and cq(0, 0^). 



7 



3 Two-dimensional finite membrane 



3.1 Coefficients Aey. 

To calculate Acy we should define explicitly what A<E>( y ) is. We recall that 
these vectors were introduced in order to write correctly the expression (^) for 
the membrane's points. In other words, relations (^) are satisfied automatically 
for any external point, but they should be imposed artificially for any point of 
the membrane. 

For membrane's point we have 

0, if (x,y)eM\S, (13) 

because particles cannot penetrate in the depth of the membrane. The first line 
of (^) represents the condition Px m .y m { n ) = &m,n- Indeed, if a particle starts 
from the boundary point (x m ,y m ), it should be immediately absorbed by this 
point. In other words, the probability to be absorbed by the m-th site is equal 
to 1 while the other sites have no chance. Expression ( |l3| ) means that in A& y ) 
there is only the contribution of boundary points and of membrane's points near 
boundary. A simple verification shows that 

(A& v ">) x , = XM(x',y') Ae lJ{x '< v ' )9 - e iJ ( x ' +1 ' y 'V (14) 

_ e iJ(x>-l,y')6 _ e iJ(x',y'+l)0 _ e iJ(x>,y'-l)6 _ ^,^(0) 

where we enumerate near-boundary functions with the help of function J(x' , y') 
(we use notation <j)[ m ] for the near-boundary function with index to). Again we 
insist to use our convention about S(x',y'), i.e., in expression @J) there are 
only the terms where corresponding point lies on S. 

Using (|l4"|), we determine Acy according to its definition, 

L 

x' — — L 

Usually there are a few nonzero components of A$' a ' for each y' . Indeed, 
according to the formula (^J), (A& v ^) x > is defined by the boundary points and 
by the membrane's points near the boundary. But on the n-th level there are a 
few such points, if n =/= 0. On the contrary, on the level zero there is an infinity 
of the boundary points due to the plane "tails". Thus, the vector (A*!^ 1 )) has 
exceptional structure. It contains the usual terms due to the non-trivial part of 
the membrane, and the contribution of plane "tails" . 



3.2 Boundary points' contribution 

Let us calculate the contribution of the m-th boundary point to 4> x ,y{6) (see 
the second term in (|12|)). During these calculations we consider the case y > 0, 
i.e., we write just omitting 7_^/ • The opposite case y < is obtained by 
reflection of all ordinates with respect to the horizontal axis, y — > —y, y' — > —y'. 
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Figure 4: Boundary points' contribution. The boundary point B = (x m ,y m ) 
with the outer normal (6x m , 8y m ) = (0, 1) is arounded by two boundary points 
on the same level {B\ and -B2), and by one external point E and one internal 
points I. 



Following the direction of the outer normal, there are four possible positions 
of the boundary point. Let us consider the case when (Sx m , Sy m ) — (0, 1) (see 
Fig. @). We suppose that the boundary point does not lie near the corners (i.e., 
we consider general position of the boundary point). The points near corners 
will be considered separately in Appendices |9.l] . Note that these points give the 
additional corrections due to the lattice distortion caused by the procedure of 
corner's removing. Usually their influence is negligible, and we do not account 
them here (see Appendices for more details). 

The boundary point B = (x m ,y m ) has four neighbours: one external point 
E, one internal point /, and two boundary points B\ and Bi with indices m\ 
and 7712. Let us write accurately the corrections for all these points. Using (|l4|), 
we obtain 



e-*(*«.-i)0» 7 ^(^). 
e -i( Xm +i)e h7 (y)( &h y 

e- fam9 *7£!-i(0h). 

The third column represents the weight factor of corresponding point (see for- 
mula (|l2|)). Now we group together all terms containing factor e™ e , and call 
such a group the boundary point's contribution Z m (we also add to this group 
the corresponding near-boundary function 4>[m]), 



B : 




= 4e tme - q 


\ m] {6) - e lmi6 - e lm ^ e 


B x : 




— 4gimi0 _ 


M ( Q\ ySmO Jim 6 
0[mi](0J - e - e 


B 2 : 




_ 4gim 2 e _ 


✓A ( (X\ siimO „im" 
0[ma] W ~ e ~ e 


I : 


A^., fc -i(J) 






E : 




= 0, 





X < 



(4-2cos0 h )7 : 

(4-2cos0 h )7 : 
( 4 _ e -ie h)7 (ri 

l(4-e 



if (5x m ,Sy m ) = (0,1), 
if (Sx m ,Sy m ) = (0,-1), 



~ ieh )ly2 



iv) Ay) 

<Vm— 1' 

• u) _ 7 (w) ; 

7^J-i " T^+i, if (fe m , ^ m ) = (-1, 0) 
7^-1 " 7^+i, if (fe m , Sy m ) = (1, 0) 
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(the first line corresponds to Fig. || three other cases can be obtained by anal- 
ogy). We see that these contributions are similar. We complete each one by 
substracting and adding corresponding term, i.e., 

e rm0- tXm e hl (y) + ^ if (fam; Sym) = (0; 1); 

g _ ^(complete) _|_ J e*™^ ix mSh _ 1 , if (5x m , 5y m ) = (0,-1), (15) 
e ime -"™^7^e^, if {6x m ,6y m ) = (-1,0), 
e imfl-ix ra fl h7 (») e -tfl hi if (Sx m ,Sy m ) = (1,0), 



(4-2co S ^) 7 to)~ 7 ^_ 1 - 7 ^ ) +1 - e -^%(^ [m] (0). 



where 

rf(complete) im6 — ix m dh 

Using the second identity of (|llj), we write 

ryl complete) imB — iXmOf, r —ix^Qy, (y) i { n\ 

z ,n -e "Oj/,j/ m -e m h 7^</> [m ](6l). 

With the help of the outer normal {8x m ,8y m ) 1 we can simplify (p"5|), 

(16) 

3.3 Distribution of hitting probabilities 



Now we can come back to the main problem. Expression (12) contains two 

terms which can be denoted as 4>x,y{6) and (j>x}y{9). We are going to simplify 
each of them independently 

Let us consider the first term, 



&0)= E ^7— tV IwI (WM/,)- 

h=-L 

Following the definition of c y (9, 9h), we have 

L 

c (9,e h ) = ($ ,v h )= e- lx ' e ^ X 'A0)- 

x' — — L 

The level zero can contain external and membrane's points. Characteristic 
functions 4'x,o{9) for the membrane's points are given by (|l|). The external 
points on the level zero were called ground points, and they are enumerated by 
index g — 1..G. Thus, 

L G 
x' = -L g=l 

Substituting this expression into </4*y(#) and taking the limit L — > oo, we obtain 

oo G 

C(«)= E ^^+£^,^,,0(9). (17) 

x' — — oo g— 1 
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This is the contribution of the level zero. 
The second term, 



n 



h=-L 



y' = -N* x' = -L 



represents other levels except level zero. Summation over x' and y' can be 
replaced by the sum of the boundary points' contributions, 



A I 



h=—L ' m=l 

Using (^6|), we write the limit L — > oo for the case y > 0, 
M 



d(9' 
2^ 



E< 

m— 1 



6 Av) 



Vm+Sy, 



l\ p ~iSx m 9' _ (y)/nl 



(0')e 



Changing the order of summation and integration, we obtain 

M 



e °x,x m o v ,y m +e u x _ Xm _ Sxm u x _ Xm <p[ m ](v ) 



(18) 



where we introduced coefficients 



(10 



(-v), 



2tt 



Note that function 71~/ (0) is added into the definition of D^' v in order to 
generalize it for the case y < 0. It means that formula ( fL8| ) is already valid for 
the both cases y > and y < 0. Using definitions (|^) and (|l0|), we rewrite this 
expression in terms of H^, 



D v ' v = 



mia{\y\,\y'\} . , 
3=1 

o, if y ■ y' < o. 



(19) 



We bring expressions ( p~7| ) and ( |l8| ) together to obtain 

M G 



where 



OG 



9=1 



a _i_ r> v > Vm+Svm 

u y,y m t ^ x - x s x „ 



(20) 



(21) 



m— — oo 



Here we should explain how ( |2l| ) is obtained. This sum over m can be 
separated in two parts: the sum over m € [1, M] and the sum over m ^ [1, M]. 
The first part contained in (pq) corresponds to the boundary points lying on all 
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levels except level zero. The second part gives contributions of boundary points 
on the level zero, including the plane "tails". In the case y > 0, expression ( p"7| ) 
contains the sum 

oo 
x' — — oo 

which can be represented as 



e imf> H. 



For such points y m = 0, Sx m — 0. Supposing 8y m — 1 and using relation 
D^' 1 = H y (for y > 0), we immediatly obtain 

m0[l,M] 

that proves representation ( pl| ) for 4> x ,y{Q)- Now, if t/ < 0, we should not write 
the contribution of plane "tails" because it is compensated by Ac_i (see the 



end of Section 3.1). In other words, the plane "tails" have no direct influence on 
the points in the lower half plane^J. The same concerns the case with y > and 
5y m < 0. In the case y = coefficient H y becomes 5-symbol which appears in 
expression ([H]) explicitly. We conclude that expressions (|2(]) and ((2^) are valid 
for any external point (x,y), in particular, with y < (if such points exist). 

Applying the inverse Fourier transform to (^0|), we obtain the distribution 
of hitting probabilities, 

M G 

Px, v (n) = 6 x , x J y , yn +Dl^f Xn - £ D^ m P [m] (n) + ^2H^ g P Sg>0 (n). 

(22) 



3.4 Equations for near-boundary and ground functions 

To complete our calculations, we should find the near-boundary functions P\ m ] ( n ) 
and ground functions Px g ,o{ n ) entering in expression (p2|). After that, one can 
use this expression for any x, y and n. We take 

x = Xk + Sxk, y = yk + Sy k for k€[l,M] 
to obtain M linear equations for the near-boundary functions Pru (n), k € [1, M] 

M G 
f[ k ][n) - ^ Xk+ Sx k -x n ~Sx n U x k +8x k ~x m ^[rn]\ n )+ 2^ n x k +Sx k -x g ^ a fi{n). 

m=l 3—1 

(23) 

Let us introduce matrices, 

(D NN ) k , m = DllXtl'^ k e [1,M], m e [1,M] 
(DN G )k,m = KlX tZ-*^ ke[l,M],me[l,G}. 

5 See Conclusions for more detailed discussion on this topic. 
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If there is no ground functions, i.e., the non-trivial part of the membrane com- 
pletely lies in the upper half plane, we take D^ G = 0. 
Now we can rewrite ( p3| ) as 

(I + D NN ) ... -D NG ... = 

V P m (n) J V P SG .o(n) J V Dll\tV-t- V L n 

(24) 

If ground points exist, we take the discrete Laplacian equations (||) for these 
points in order to obtain conditions for ground functions, 

iP Xg ,o(n) - Px 9 -i,o(n) - Pz 8 +i,o0) - P* B ,i(n) - Pa„_i(n) = 0, g£ [1..G]. 
Substituting Pg s ,i(n) and P Xg ^i(n) from ( p2| ) into these conditions, we obtain 

M G 

- E (^f-x m + 4;-l m m )PM(n) + £ 2//,'.. ,../',.,..,:«:. 

771=1 ff 7 — 1 

In matrix form, 

^u(n) \ / ^i,o(») \ i / ^l;i 5 !- fa „ 

#Giv| ... + #GG =5 ... |, (25) 

P[M)(n) J V P *cAn) J V fl^tt*" 



— a:^ — (5x r , 



where 



{2, if fc = to, 

-0.5, if |x fc -S m | = 1, k£ [1,G], me [1,G], 

0, otherwise 

(Aav)*,*, = ^I:U m , * G [1, G], m e [1, M]. 

We have two matrix equations ( |2~i| ) and ( p5| ) which allow to find near- 
boundary and ground functions, 

P[i](n) \ ( Pfi](») 

... =(/+r ... ! (26) 

P[M](n) J V P [V]M 



where 



£> = + D NG D GG D GN and 



-P* (Til \ / Tjy^+^y^'Vn+^Vn \ / TTVn+Syn 



+ -D NG D GG 

P* (n) J \ nVM+8vM>vn+8u n J \ Tjy™+ S yi 



(27) 



Using ( p5|) again, we obtain ground functions, 

ry n +$y 



= ^gg - - J D G ^ GJV (/+^)- 1 

(28) 



13 



Expression ( p2| ) is our main result. What have we done? Using the char- 
acteristic functions technique, we express the hitting probability P Xjy (n) (for 
any x, y and n) in terms of the explicit coefficients D^' v and finite number of 
coefficients P[ m ](n) and P Xg ,o(n) which can be calculated with the help of (p6|), 
© and (||). 

Note that P x , y (n) depends on n only through coordinates x n and y n of the 
n-th boundary point and through outer normal (5x n , 5y n ) at this point. It 
means that obtained distribution of hitting probabilities does not depend on a 
choice of parametrization of the membrane. In other words, we can use any 
parametrization of the membrane's boundary. This remark will be used in the 
next section. 



4 Multi-dimensional membranes 

We confined ourselves to the two-dimensional case in order to clarify the treat- 
ment. Now we are going to generalize the previous results for multi-dimensional 
membranes. 

Consider d-dimensional hypercubic lattice. As above, we can define a finite 
membrane with the help of similar conditions as in Section |j In particular, 
we suppose that membranes have a compact support (the compactness condi- 
tion), i.e., there exists a (d — l)-dimcnsional hyperplanc (level zero) such that 
the membrane's boundary is just its "finite" perturbation. We choose such co- 
ordinates that this hyperplane is defined by equation x^ d > — 0, its points are 
enumerated by multi-index x = {x^ x \ x^ -1 '). Coordinate x^ d ' of the orthog- 
onal direction is denoted y. Let us denote e& the unit vectors in k-th direction, 
k S {l,...,d- 1}. 

As usual, we enumerate the boundary points by index n, (x„, y n ). The outer 
normal on n-th point is (Jx„, £?/„). We can introduce the n-th level as a set of 
points lying on the hyperplane with y = n, i.e., {(x, y) : y — n, Vx}. 

The discrete Laplacian equation in d-dimensional case is 

Px,y( n ) = ^5Z( Px+e -^( n ) +- p x-e 4 , 1/ (n)J + — (Px, y+ i(n) + P*,y-i{n) 

whence one can obtain the distribution of hitting probabilities for a planar 
multi-dimensional case, 

Hl= J ...J j^^ e <C**y»l( g< ) (29) 

— TT — TT 

where (x • 9) is a scalar product of two (d — l)-dimensional vectors x and 8; 
function tp(6i) depends on 0\,..., 6 c i-i, 



-l 

cos( 

1=1 



\ 



d-l s 2 

d — cos 8i J — 1 



The crucial idea of our generalization is that all manipulations along the 
orthogonal direction remain valid. In particular, the expressions (^) and ( [l0| ) 
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are true, if one takes f{6\, 6d~i) instead of (p(9). It means that we can impose 
the same conditions (||) basing again on the technique proposed in . Now one 
can repeat the previous calculations of Section |3|: 

- to calculate coefficients Ac y ; 

- to obtain contributions of boundary points; 

- to impose conditions for near-boundary and ground functions. 

Having made these technical steps, we understand that the only coefficients 
H.% (and D% v as consequence) are different for d > 2, but the structure of 
solution is exactly the same. It means that we leave expression ( p2|) without 
changes for d-dimensional case, 

M G 
m— 1 y— 1 

(30) 



Here coefficients -D^-x are s * n ^ defined with the help of (19) but for coefficients 
Hx-x' we should use ( p9| ) instead of (||). 

As above, we can obtain a set of linear equations for near-boundary and 
ground functions. We take 

M G 

^[fe](n) - ^x„+Sx k -x n -Sx n - 2^ ^x fc +5x fc -x m ^[m](") + Z^^x fc +A-x fc -x g ^x,,0(n) 

m=l (?=1 

for near-boundary functions, and 

1 d_1 1 
^,oW = —J2(Px g+ eM + P* g - ei , (n)) + 2j(Ps B ,i(n)+^*.,-i(»)) 

i=l 



for ground functions. Using the same representation as in Section 3.4, we obtain 
the near-boundary and ground functions according to expressions fl26|) and (p8|). 
The only change that we should make is concerned with matrix Dqq. Indeed, 
we generalize this matrix as 

{d, if k = to, 
-0.5, if |x fe -x m | = 1, 
0, otherwise. 

It is important to stress that we used the remark at the end of Section [| the 
distribution of hitting probabilities should not depend on a choice of parametriza- 
tion. It allows to avoid a complex multi-dimensional paramctrization of the 
membrane. To be more rigorous, one should introduce such a parametrization, 
recalculate distribution P x<y again, and re-enumerate the boundary points by 
the single index in order to obtain (|3C|). We omit these technical details. 

Note that expression (^0|) has the same structure for any dimension of the 
lattice, the only distinction is contained in coefficients HQ which are individual 
for each d. We give a useful asymptotics for coefficients HQ, 

Y{d/2) \y\ 

n d/2 ( x 2 +2/ 2)d/2 V ) 

that is the multi-dimensional generalization of the well-known Cauchy distribu- 
tion of hitting probabilities of the Brownian motion. 
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5 Brownian self-transport operator 



As it was mentioned before, we are interested in the problems of Laplacian 
transfer, and consequently, in the Brownian self-transport operator Q. We recall 
that Qk, n is the probability that a random walker contacts at the first hit the 
n-th site of the membrane if started from the fc-th site, without touching the 
other sites. Obviously, any path starting from fc-th site of membrane must pass 
through the corresponding near-boundary point. Provided that we removed all 
the singular (corner) points, any boundary point has exactly one near-boundary 
point. It means that 

Qk,n = P[k](n), 

i.e., the matrix Q is composed of the near-boundary functions P\j^{n). If one 
denotes 

Qk,n = -P[fc]( n )i 

expression (^6|) becomes 

Q = {I + D)- 1 Q*. (32) 

Recall that the system of linear equations ( p3| ) was obtained only for non-trivial 
part of the membrane, in particular, matrix Q has MxM elements. However, we 
are usually interested in the whole Brownian self-transport operator, including 
the plane "tails". 

To construct the Brownian self-transport operator for all sites of our interest, 
we can slightly modify definitions of matrices Dnn, -Djvg and Dqn- Let A be 
the set of indices of sites for which one wants define Q. For example, if one 
needs to know Q for all sites of the membrane, it is sufficient to take A = Z. 
For the numerical simulations one takes a finite number of sites. In any case, A 
must contain all non-trivial sites, i.e., [1..M] C A. We remark that summation 
over m in ( p2| ) can be elarged to all integer numbers from -ootooo. Indeed, for 
m £ [1,M] we have y m = 0, and Bjj.'*™ = according to (|l^). In particular, 
this summation covers all possible sites which can be of our interest. Now we 
rewrite the definitions of matrices Dnn, Dng and Dgn, 

k e A, m e A 
he A, me [1,G\. 
k G [1, G], me A. 



N N )k,m - U Xk+Sxk _ Xrr 
(D NG )k,m = H%" h Xsl k k -x„ 

{DcN)k, m = \HtT-x m : 



Matrix Dqq remains unchanged, and 

D = Dnn + DngDgg^gn- 

Evidently, matrix Q* (i.e., vectors P\ m y) should be also recalculated for all sites 
of the interest. 

These modifications were made to obtain the Brownian self-transport oper- 
ator in the matrix form ([52]). We stress that after finding the near-boundary 
and ground functions, one can use ( |22] ) directly to "compose" the matrix Q, 

M G 

^k,n - LJ x k +Sx k -x n ~Sx k U x k +&x k -x m ^{r<A\n) j- } y -"g fe -x n ^x g ,o(n). 

m—1 g—1 



1G 



We can conclude that the Brownian self-transport operator is explicitly con- 
structed for any reasonable multi-dimensional finite membrane. Now it can be 
useful to study the general properties of Q* and the spectral properties of D, 
and their dependencies on the geometry of the membrane. This work is not yet 
finished. 



6 Generalizations 

As we explained in Section H, expression ( p2| ) can be generalized for the multi- 
dimensional case. One can go further by regarding more complex problems. 
In this section we consider the influence of different barriers, and the time- 
dependent distributions of hitting probabilities. 

It is important to stress that we have obtained the Brownian self-transport 
operator for a membrane with infinite "tails" placed in the infinite space. How- 
ever, usual physical membranes have a finite size, and they are placed in a closed 
volume (like a box or a cell). In order to apply our results to real membranes, 
these features should be taken into account, for example, by introducing absorb- 
ing or reflecting barriers. First, we consider the horizontal barriers, and then 
the vertical barriers. 



6.1 Horizontal barriers 

Suppose that there exists a horizontal barrier at the level JV + 1 > N which 
absorbs or reflects all particles. In this case, it is not difficult to recalculate 
the distribution of hitting probabilities (|2^). Let us return to the recurrence 
relations (Q). In order to close them, one can use a condition like 

with a certain coefficient 77 which can depend on 9 and Oh- Now one can express 
c v in terms of Co and {Ac,/} (and of some explicit functions like 77 and if). This 
expression was found in |7j] , 

('11 



Cy = .f^(c +J2a y/ Acy / )+ay £ /^Acy, 
^ y'=l ' y'=y+l 



where 



and 



fW = [1 - wm ^-y^{6 h )M0 h ) - rj\ 

J y [1 - W (e h )} - <p^+i{o h )\<p{e h ) - v ] v K h> 

When level M of the barrier goes to infinity, function fy^ tends to fy°°^ = ip y , 
independently of the value of rj. It simply means that one can choose any 
condition at infinity. 
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The absorbing barrier corresponds to r\ — since the particles cannot reach 
the membrane from any point (x,J\f + 1); thus 

i , n 2SS+2-2y 

The reflecting barrier corresponds to n = 1 since the particles should have the 
same probability to reach the membrane from points (x,Af+l) and (x,Af); thus 

i I ,JW+l-2y 

Jv l + (/3 2A^+l V ■ 

All other calculations are left unchanged, but we should replace ip y by f^- abs ^ 
or ff< ref) . It means that one can use expression (E2h where coefficients H% are 
replaced by modified ones, 

f if y > TV, 

H y x = I J f> e '#V) if 1 < V < JV, 

I " if y < 0. 

The last line represents the fact that the barrier in the upper half plane has no 
influence on the lower half plane. 

Using this definition of coefficients H% , one can calculate the distribution of 
hitting probabilities for a general membrane with an additional barrier absorbing 
or reflecting particles. Note that these modifications are valid for any dimension. 



6.2 Vertical barriers 

Now we can limit the particle's movement along the horizontal axis putting an 
absorbing or reflecting vertical barrier. Let us put two barriers at x = — L — 1 
and x = L + 1 so that the non-trivial part of the membrane lies between these 
vertical lines. In Section ^, we represented the Laplacian equations (^) in matrix 
form using matrix A, 

A M =4, A i>i+ i = A i+ i,i = -1, A- LtL = A L - L = -1, 

and all other elements are 0. The last equality was imposed artificially to 
obtain the cyclic structure of matrix A that was convenient to have explicit 
eigenvectors Vh- After all calculations, we took the limit L — > oo, there- 
fore this modification vanished. Now we are working with finite L, conse- 
quently, we should be more accurate with these conditions. Note that relation 
A-lx = Al_-l = — 1 corresponds to periodic or cyclic boundary condition, i.e., 
we maintain 4>L+i,y = 4>-L,y and 4>-L-i,y = 4>L.y In other words, we suppose 
that each pair of points (L + 1, y) and (— L, y) is identified or "glued together" 
(the same for points (— L — l,y) and (L,y)). 

In order to introduce vertical barriers, one should change the boundary con- 
ditions. For two absorbing barriers one should simply remove the artificial 
condition, i.e., A-l,l = Al,-l = given that (j)L+i,y = 4>-L-i,y = (particles 
are absorbed by barriers). For two reflecting barriers one should also introduce 
the following modifications, A-l,-l = = 3 given that <j>L+i, y — 4>L, y and 
4>—L—i,y = 4>—L,y Obviously, one can also consider the mixed case with one 
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absorbing barrier and one reflecting barrier. Eigenvalues and eigenvectors of 
these modified matrices have no explicit form. It is more convenient to use the 
initial matrix A (cyclic boundary conditions) with additional corrections, i.e., 

where we introduced a new correction vector A$W whose components are equal 
to except 

(A$<»))_l = -H, v (A$(»)) L = -</>- L , v (absorption) 

(A$fe))_ L = -</> LiV + 4>- L , y (A¥^) L = -(j)-L,y + H,y (reflection). 

Now we can repeat the previous calculations to obtain a new distribution of 
hitting probabilities P x y (n) instead of (j22]), 

M G 

P x , y (n) = D™_l + n 5 _f Xn - Dl'^ m P [m] {n) + £ H y x _ Sg P Sgfi (n)~ (33) 

rn—1 y— 1 

oc , i 

E [D™ l PlA*) + D Z+ L P-L,y> (»)] , (absorption), 

y'=l 

OO , , 

E - Dlf L )[P L , yl {n) ~ P-L,y(n)], (reflection), 

y'=i 

where H y should be calculated as finite sum over h from —L to L, 

P ix °h O-Trh 



2L + l r v " 2L+ 1 



and is expressed according to (|19|) in terms of H%. 

We should make several remarks. First, for sufficiently large L one can use 
the integral expression (|^) for coefficients H y as an approximation. Unfortu- 
nately, such approximation cannot give an accurate result for \x\ near L. When 
we are not interested in points near barriers, it does not lead to the problems. 
But in some cases it can do. Therefore one should be careful working with finite 
membranes. 

Second, for absorbing and reflecting barriers there appear new unknown 
functions P±L.y{n) which can be called barrier functions. It means that one 
should write additional conditions to obtain a closed system of linear equations 
for near-boundary, ground and barrier functions. In order to obtain the con- 
dition corres pon ding to the barrier functions 4>±L,y, one can maintain simply 
x = ±L in (]33j) taking different y. The problem is that there is an infinite 
number of barrier functions. A convenient solution is to consider two vertical 
barriers together with a horizontal barrier (that is the usual case for a real 
physical membrane). In this case, the number of barrier functions is just Af. 
However, one should not forget to use coefficients H y and D y ' v corrected by 
introducing the horizontal barrier, see the previous subsection. 

Third, one can try to use a rough approximation for a very large L. Indeed, 
if \x\ -C L, one has ~ l/L 2 , and the barrier functions can be neglected. It 

simply means that one uses the initial formula ( p2| ) to approximate the solution 
with two remote barriers. 
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Fourth, in the case of one horizontal and two vertical barriers there is a finite 
number of external points, i.e., one can solve the system of discrete Laplacian 
equations (|^) properly, without using probabilistic techniques. In other words, 
the essential advantage of the present formalism is the possibility to analyse the 
membranes with an infinite number of sites. 



6.3 Time-dependent distribution of hitting probabilities 

Previous results were related to the Laplacian equation, i.e., we considered dis- 
tributions of hitting probabilities independent of time. Generalizing the prob- 
lem, we can find the distribution of probabilities Px} y {n) to hit the n-th site 
of the membrane at i-th step (here t is discrete time; for the Brownian motion 
it is continuous time) starting from an external point (x,y). From the techni- 
cal point of view, it is convenient to consider the Laplacian transform of this 
distribution, 

oo 

^(n) = EA*pW(n). 
t=o 

Note that A = 1 corresponds to the previous time-independent distribution. 
More generally, Pj A J (n) gives the distribution of hitting probabilities if particles 
have killing rate (1 — A). For the planar case, one can apply the same technique 
as usual to obtain the distribution 



d6 
— ( 
2tt 



t i(x — n)t 



^(0;A)=^_„(A) 



with function 



<p(0; A) = 2/A - cos 9 - ^(2/ X- cost 



This expression can be easily rewritten for the multi-dimensional case. Now we 
repeat the previous calculations for general membranes to obtain exactly the 
same formula as (22), 



P^H=S x , x J y , yn +D 



y,y n +Sy n 
.—5x- 



.(a)- 



M 

V D v ' Vm 

/ j X — Xn 

rn — 1 



G 

E 

9=1 



where coefficients D}j,' v (A) are computed according to formula ( fl9| ) with H}j.(\) 
instead of H%. Imposing usual conditions for near-boundary and ground func- 
tions, one can close these relations in order to find P,^|(n) and P^ <0 (ri). It is 
important that A appears in all expressions as a parameter. 
Now distribution Pi*y(n) can be found as 

1 

' dA* 



Replacing function ip(9; A) by 



t\ 



A=0 



p(0i,...,0 < ,_i;A) 




\ 



d-l 

E 

i=l 



cos 9i 



1, 



we easily generalize this time-dependent distribution of hitting probabilities to 
the multi-dimensional case. 



20 



162 9f\ 



1214 



31 



155 



3n 



62k. 



55 



93 



Figure 5: Two simple membranes : (a) convex membrane (N = 31, A^* = 0, 
M = 155); (b) concave membrane (N = 0, N* = 31, M = 155). 



7 Numerical verifications 

In this section we briefly present some numerical results to check the validity of 
the method. 

For a particular non-trivial membrane, the validity of expressions (p2| ) and 
( j30|) can be verified by comparing them with the results obtained by numerical 
simulations (simple random walks on the lattice). We have taken two- and 
three-dimensional membranes represented on Fig. ^| and Fig. [| 

In the first case (Fig. [s]a), there are no ground functions, G = 0, and there are 
155 near-boundary functions which can be calculated using (j32|). For the con- 
cave membrane (Fig. |]b), there is 31 ground functions and 155 near-boundary 
functions. On Fig. (?]we present the distribution of hitting probabilities P{2o\(n) 
(i.e., Q20,n) obtained with the help of formula (|2^). The solid line corresponds to 
the convex membrane, the dashed line - to the concave membrane. These results 
are obtained for the membranes placed in the finite box. In other words, we take 
the horizontal absorbing barrier on the distance D = 100 from the ground level, 
and two vertical barriers under the cyclic boundary conditions (see Section ^|) . 



Moreover, we account the corner points' corrections (see Appendices 9.1). 

Fig. U shows the distribution of hitting probabilities P[ix3](n) (he., Qn3. n ) 
for three-dimensional membranes (see Fig. ^|): the solid line corresponds to the 
convex membrane, and the dashed line - to the concave membrane. These 
curves are obtained with the help of formula ([52]) , where we used (|2^) for H~[ 
instead of @. We calculate these distributions for the membranes placed in the 
finite box which is formed by the horizontal absorbing barrier at D = 15, and 
by four vertical barriers with cyclic boundary conditions (see Section |^) . Also 
we took into account the corner points' corrections. 

In order to confirm that the method works, we performed the Monte-Carlo 
simulations. Starting from the near-boundary point {xk + yu + <%k), the 
particle walks at random on the lattice until it is absorbed either on the mem- 
brane's boundary or on the horizontal barrier. Repeating this procedure, we 
obtain the numerical distribution of frequencies of visits for all boundary points. 
Dividing by the total number of walks A^, we calculate the approximate hit- 
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Figure 6: Three-dimensional membranes: (a) convex membrane (N = 6, N* = 
0, M = 365); (b) concave membrane (N = 0, TV* = 6, M = 365). There are 
13 facets 5x5 (including 8 facets on the ground level), and 8 segments 5x1. 
Boundary points are enumerated by index n following the rules: 1 st facet (y = 0): 
n = (\z\ - 1) -5 + x, 2 nd facet (x = 6): n = (\z\ - 1) ■ 5 + y + 25, 3 td facet (y = 6): 
n = (\z\ - 1) ■ 5 + x + 50, 4 th facet (x = 0): n = (|z| - 1) • 5 + y + 75, 5 th facet 
(\z\ =6): n = (y -l)-5 + x+ 100, 6 th - 13 th facets (z = 0) - by the same way. 




Figure 7: Distributions of hitting probabilities (in log-scale) Qk,n with k = 20: 
the probability of the first contact with points (x n ,y n ) of membrane's boundary 
if started from the point (x 2 o, 2/2o)- 
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Figure 8: Distribution of hitting probabilities (in log-scale) Qk, n with k = 113 
(the center of the fifth facet) for the three-dimensional convex and concave 
membranes. 

ting probability Qk.n for each pair of sites k and n. When Nt goes to infinity, 
the approximate values Qk, n tend to the exact values Qk, n - These simulations 
represent the well-known Bernoulli trials. In order to estimate the accuracy of 
values Qk.n for the finite Nt, we use the central limit theorem. It gives the 
natural measure AQk, n of deviations Qk.n — Qk,n, 

AQ k ,n ~ y/Qk^/y/Ni, 

where a ~ y/Qk,n i s the dispersion of the Bernoulli trials. 

Taking different values of Nt and performing the Monte-Carlo simulations for 
all possible sites k and n, we compute the maximal deviation of the approximate 
value Qk, n from the exact value Qk,n, 

{\Qk,n-Qk,n\ /rrl 
— = J — V N t > ■ 
\/Qk,n J 

Values of a for different membranes are given in Table 1. 



Membranes 


10 5 


10 B 


10 v 


2D convex 


4.99 


4.02 


4.43 


2D concave 


4.96 


4.91 


4.46 


3D convex 


5.33 


4.44 


4.46 


3D concave 


6.33 


4.59 


4.46 



Table 1. Monte-Carlo verifications of the method's work. For N t £ {10 5 , 10 6 , 10 7 }, 
the maximal deviation a is calculated. 

These results confirm that formulae ( p2|) and (|30| ) work correctly. 
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8 Discussion and conclusions 



In this paper, we defined a wide class of finite membranes which play an impor- 
tant role for the discrete Laplacian transfer problem. We based on the technique 
of characteristic functions developped in our previous paper. In particular, ex- 
pressions (||) are used to close the recurrence relations (0). We remind that 
expression c/v„+i ~ yc/v^ was the central approximation of H. In order to 
obtain this approximate relation, we supposed that from a remote line y = N u 
(for large N u ) the membrane can be viewed as almost translationally invariant 
object (along the horizontal axis). Then the explicit solution of recurrence re- 
lations (Q) had been found. Finally, we took the limit N u — ► oo (and Ni — > oo 
for the lower half plane) to obtain expression (w for coefficients c y (Q,9h). It 
means that influence of approximate relation (0) vanishes. In other words, the 
approximate relation (^) is taken at infinity (N u — » oo), i.e., this condition has 
no influence on the solution. Moreover, one can verify directly, without any 
approximation, that (JsJ) is the exact solution of recurrence relations (^) . 

We stress that the essential result of our calculations is the form of solu- 
tion ( p^ ) which expresses the distribution of hitting probabilities for a general 
membrane, P X y(n), in terms of corresponding planar distribution, H%. In other 
words, we explicitate how the membrane's geometry should change the planar 
distribution. It means that having solved the planar problem, one can easily 
generalize its solution for the case with a rather complex geometry. For exam- 
ple, when the problem with horizontal barrier for the planar case was solved, we 
just replaced old coefficients H% by the new ones. The same ideas were used for 
the time-dependent distribution and for the multi-dimensional generalization. 
Note that these motivations are frequently used in the theory of conformal 
transforms. Indeed, taking an appropriate conformal transform, one can map 
the initial complex set into a simple set (like a half plane or a disk), solve the 
simplified problem, and then reconstruct the solution by the inverse conformal 
transform. Note also that the conformal theory works only in two-dimensional 
case while the present approach is valid for any dimension of the lattice. 



8.1 Ground functions 

Let us briefly discuss the role of ground functions. From the beginning, the 
upper and the lower half planes were considered separately (for example, see two 
terms 7^, and 7!^/ in equation (||)). The technique proposed in (7) consists 
to step down from the N u -th and (— N{)-th levels to the level zero. It allows to 
express coefficients c y in terms of Ac y . But we cannot step down directly from 
7V„-th to (— iVj)-th level, because there appears an infinity of near-boundary 
functions due to (A<E>o), and we do not achieve our goa|| In other words, we 
cannot "pass" through the level zero in our treatment. On the other hand, a 
random walk started at the upper half plane may hit the boundary point in the 
lower half plane (if such point exists). This walk passes the level zero through an 
external point, i.e., through a ground point. Therefore we can say that ground 
functions "connect" the solutions for the upper and lower half planes. 

Consider as example the concave membrane (Fig. ||d). The probability to 
hit the n-th site in the lower half plane if started from a point (x, y) in the upper 

6 For the same reason the compactness condition was imposed at the beginning. 
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half plane becomes 

G 

9=1 

We simply took ( p2| ) under condition that y > 0, y m < for m <E [1, M]. This 
formula has a clear probabilistic sense: to hit the n-th site, the walk should 
reach one of the ground points (enumerated by g) with probability^] H%_ i , and 
after that hit the n-th site with probability Px g ,o(n). 

8.2 Probabilistic sense of coefficients D% v 

Now we clarify the nature of coefficients DV' y which seem to be artificial at the 
first sight. First of all, they have no purely probabilistic sense (for example, they 
can have values greater than 1). To understand the sense of these coefficients, 
consider the simple problem: to find the probability to hit a point (xo,yo) if 
started from [x, y) without touching the horizontal axis. The solution of this 
problem is given in Appendices 9.2, 

P[(x,y)-*(x a ,y )] = n ya,yo ' (^) 

Using this solution, we can rewrite ( p^ ) for (x, y) ^ S as 

Px,y(n) — P[(x,y)~*(x n +8x n ,y n +8y n )]DQ n+5y "' yn+Sy " (35) 
M G 

~ H P l(x,y)^(x m ,y m )] D o m ' Vmp [m](n) +J2H y ^ Sg P Xgfl (n). 

Now this expression has more clear sense: we take all possible ways to reach 
near-boundary point (x n + Sx n , y n + 6y n ) without touching the horizontal axis, 
and then subtract all ways which go through membrane's points. Coefficients 
jjy m ,y m anc j £)y n +Sy n ,y n +Sy n kave combinatoric origin: they allow to account 
different possibilities to go through the membrane. According to definition 
©, 

J'=l 

and for two-dimensional case the approximation 

D y ' v « — — + const (const w 0.3675) 

gives Dq' v with a sufficiently high accuracy. 

Writing ( |35| ) for the planar membrane, one has 

Px, v (n) = P[(x,y)^(x n ,i)]Dl' 1 , 

i.e., aparticle reaches near-boundary point (x nj 1) with probability P\( x ,y)^(x n ,i)]i 
and then it hits corresponding boundary point (x n ,0) with probability Dq' 1 = 
Hq just as required. 



7 We write 

membrane completely lies in the lower half plane. 
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8.3 Coefficients D™ x , like Green's functions 

Coefficients D x ^ x , can be viewed from another standpoint. A direct verification 
shows that 

L ' 1J x-x l — u x,x'Vy,y', ^ x-x' — u ' 

where the discrete Laplacian L operates on coordinates (x,y) or (x',y') (we 
remind that t = D v x ,^ x ). For example, in two-dimensional case we have 

^tlx^y — 'lu^y Ux-\-X } y ^x— l,y ^x.y+1 ^x,y — 1- 

It means that coefficients D x '^ x , can be treated as Green's functions of the 
discrete Dirichlct problem in the upper half plane. Consequently, the solution 
of the general Dirichlet problem in the same (planar) geometry, 

£v>x,y — fx.yi V'X.O 0, 

is given as 

u x,y = D x -x' fx'.y' ■ 

x',y' 

Unfortunately, this result has no direct connection with our actual problem 
of general membranes. Nevertheless, we can use the following trick. What 
violates the using of Green's function D x ^ x ,7 The answer is that there exist 
the boundary and internal points of the membrane with y > 0. Obviously, the 
Laplacian equation becomes invalid on these points. However, we can correct 
this situation by introducing a certain function f x ,y(n) (here n is a parameter). 
In other words, we impose the Laplacian equation artificially for the boundary 
and internal points. The correction function is equal to for any external point. 
Moreover, f XjV equals to for almost all internal point except the internal layer 
whose points are nearest to the membrane's boundary. For the boundary points 
we have 

fx miVm (n) = 4Px m ,y m (n)-P Xm -i,y m {n)-P Xm+li y m (n)-P Xmiym+ i(n)-P Xmi y m -i(n), 
whereas for the nearest internal layer 

fxm,—&X< m ,y m , — dy m \Jl) ^CmjJ/ml^)' 

Now we can write the general solution in terms of Green's functions, 

M 

Px, v (n) = E [ D x-xJ^v,An) + Dl^J^J Xm _sx m , ym -s Vm (n)] ■ 

m— 1 

Imposing the boundary condition P Xm ,y m {n) = S m , n and considering step by 
step four possible directions of the outer normal, one can demonstrate that this 
formula can be written as 

M 

Px, y (n) = 5 x , x J y ,y n + Df*£%l n - £ Dl^ m P [m] {n). 

m=l 

Note that here we do not discuss neither contribution of plane "tails" , nor 
ground functions. One can easily complete these motivations in order to obtain 
formula (^3). This approach is also valid for multi-dimensional case. 
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8.4 Comments on the Brownian self-transport operator 



The important application of the distribution (]22j) of hitting probabilities is the 
explicit construction of the Brownian self-transport operator Q which governes 



the Laplacian transfer. We obtained the explicit formula (32) for Q, and it opens 
the possibility of analytical researches in this field. Note that exact analytical 
expression of Q even for a planar membrane allowed to obtain the important 
characteristics of the Laplacian transfer, ||. 

On the other hand, expression ([32j) simplifies also the numerical treatment 
of the problem. Usually one uses computer simulations of random walks to 
calculate the matrix elements of Q. In order to obtain the whole matrix Q 
with a high accuracy, one should make enormous number of random walks. 
Moreover, one walk can be very long, especially on lattices with d > 2. On 
the contrary, working with our approach, one just needs to manipulate in the 
framework of linear algebra. Once calculated, coefficients can be easily used 
for any membrane. Therefore, one just needs to "compose" the matrix D for 
a given membrane, and to inverse (7 + D). The time required to make these 
operations depends only on the number of sites for which one calculates Q. 

In order to adjust the formalism to a more realistic membrane, one can 
introduce different barriers. As we described above, the horizontal barrier can 
be introduced by simple modification of coefficients H^, while the problem with 
vectical barriers is more complicated. 

We conclude that the present treatment opens a wide field for further inves- 
tigations. 



Acknowledgement 

The author would like to thank Professor B.Sapoval and Professor M.Filoche 
whose invaluable advice helped to develop the present approach and to construct 
the Brownian self-transport operator. 



9 Appendices 

9.1 Corner points' corrections 

In Section [3.2] , we calculated the contributions of the boundary points supposing 
that these points are in general position. In other words, we did not consider 
the neighbourhood of the corners. Moreover, the procedure of corner points' 
removing distorts the lattice near corners (see Fig. ||d), and we did not ac- 
count the influence of such distortion. For sufficiently regular membranes, these 
corrections should not change the distribution of hitting probabilities globally. 
However, the lattice distortions can perturb significantly the solution locally, in 
the vicinity of the corners. It can be important for certain problems. Here we 
briefly explain how solution ( p2|) can be improved in order to account such local 
distortions. Note that these corrections carry a particular character, they de- 
pend on the way how the lattice has been distorted. For example, if one would 
not like to remove the corners, there is no correction. 
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(x c -l,y c ) 




r J\x c ,y c -l) 



Figure 9: Corner points' corrections. 



Why problems appear near the corners? Consider as example the corner 
point (x c , y c ) shown on Fig. Ej. 

We wrote the Laplacian equation (Q) for any external point, in particular, 



0. 



After that, we removed the corner point (x c ,y c ) (thus one can take 4> Xa ,y c — 0), 
and connected points {x c ,y c + 1) and [x c — l,y c ) to each other by a liaison. It 
means that the true Laplacian equation is 



»x a ,y c + l 



9 z c + l,y c +l 



Z>x c ,y c +2 



J x c -l,y, 



= 0. 



We prefer to work with the previous general form, and for these purposes we 
introduce here the correction term A0 Xe .j, c+ i which is equal to 4> Xc ~i,y c - The 
same concerns the point (x c — 1, y c ). So, each external corner gives two correction 
terms (in two-dimensional case). Note that these corrections have not been 
accounted above, they appear due to distortion of the lattice. We should also 
consider the Laplacian equation at point (x c , y c ) because this removed point has 
been accounted in the general treatment, 



x c,yc 



Px c +1,K C 



J x c ,y c + 1 



J %c,y 



Above we supposed that the left hand side is equal to because the point 
{x c , y c ) was referred as external point. Now we calculate the corrections due 
to the lattice distortions, and for these purposes we write the correction term 
A(j) Xc} y c . Taking all these terms, we obtain the whole contribution of the corner 
point (x c ,y c ), 



Z - - \ P - ix " e h (y)_ -ix c e h Av) i ± 



L e lyc c lyc IVxcVc+l 



_JJ{x ,y c -l)e -ix c 9 h (y) _ ij(x c +l,y c )9 -ix c 9 h (y) 

Two last terms should be included into the contribution of boundary points with 
indices J{x Cl y c — 1) and J(x c + 1, y c ) respectively (exactly these two terms are 
missed if one considers the boundary points near corners, see Section 3.2). Two 



first terms involve the unknown characteristic functions <ft Xc -i.y c and <fi Xc ,y c +i 
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that can be called corner functions. It means that in general case formula ( p2| 
becomes 

corner 

f n^^ n v,vc+Sy i c 1) \ p / x 

where P xy {n) denotes the right hand side of (p2|), and (&e£ 1,2 \ fo/i 1 ' 2 ^) are two 
outer normals on the corner point (x c ,y c ). 

The corner functions P , , n,2) . r (1,21 (nl can be calculated by the same 

x a +Sx c ,yc+oy c v ' J 

way that was used for near-boundary functions, i.e., we close the system of 
linear equations and solve it. The last formula gives a right distribution for all 
sites of the membrane, including the points near corners. 



(2) 



P^+S^^+Sy^^y 



9.2 One point problem 

Consider the simple problem of finding the distribution of hitting probabilities 
Px,y{n) on the horizontal axis if there exists an absorbing point (0,yo) with 
yo > 0. As consequence, we shall find the probability to hit the point (0,j/o) if 
started from (x,y) without touching the horizontal axis. 

Here we cannot apply the technique of the planar case given that 

- there is no translational invariance ; 

- the Laplacian equation (||) is invalid on the point (0,yo); 

On the contrary, we can use formula ([|). There exists the only correction 
A$( v <>' for the unique point (0,yo), i.e., 

A Cyo(°) = -^l,yo( ) - ^-l,yo( e ) ~ ^0,2,0+1(0) - <t>0,y -l(8) 



(we omitted 40 o ,y o because it equals to 0). Using formula (|12|) and taking the 
limit L — ► 00, we obtain 

7T 

&>,y(*) = / ^e ix6 ' (y y {6')co{0,6>) + ^{e')&c yo {6) 

where 

00 

Cq(0, 6') = J2 e lx ^- 8 "> = 2ttS{9 - 9'). 



X— — OO 



Applying the inverse Fourier transform, we obtain 
P x , y (n) =Ht_ n -Dl'y°(p 1 M + P-iM + P Q ,y 0+ i{n)+P 0m ^{n) ). (36) 



Looking on this relation, we say that there are four near-boundary functions on 
four points around (0, yo). Taking corresponding x and y and summarizing four 
equations for near-boundary functions, we write the equation for their sum, 



]T P W (n) = (Hfl n + Hl\_ n + Hl^ 1 + Hl^ 
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(Df' ya + D y °i vo + D y 0+1 ' vo + D%°- hvo ) P[m](n). 



m—l 

Using the properties of coefficients H y and D y,y , we can simplify this relation 
and express the sum of near-boundary functions, 

^ P[m]( n ) = jjycyo ■ 
m=l 



Substituting this sum into (36), we have 

r>y,yo 

P (ri) = H y ~ H y ° — 

^0 

So, the problem is solved. The probability to hit the point (0,yo) without 
touching the horizontal axis is 

p(one) _ , p / s _ x 

r x,y ~ 1 r *>v\"') ~ jyyo,vo • 

n— — oo 

We hope that the consideration of this simple problem clarifies the sense of 
coefficients D y ' y and simplifies the understanding of general results. 
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